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Abstract 

The  P-vector  inverse  method  has  been  successfully  used  to  invert  the  absolute  velocity  from  hydrographic  data  for  the  extra- 
equatorial  hemispheres,  but  not  for  the  equatorial  region  since  it  is  based  on  the  geostrophic  balance.  A  smooth  interpolation 
scheme  across  the  equator  is  developed  in  this  study  to  bring  together  the  two  already  known  solutions  (P-vectors)  for  the  extra- 
equatorial  hemispheres.  This  model  contains  four  major  components:  (a)  the  P-vector  inverse  model  to  obtain  the  solutions  for  the 
two  extra-equatorial  hemispheres,  (b)  the  objective  method  to  determine  the  *F-values  at  individual  islands,  (c)  the  Poisson 
equation-solver  to  obtain  the  IT-values  over  the  equatorial  region  from  the  volume  transport  vorticity  equation,  and  (d)  the  Poisson 
equation-solver  to  obtain  the  *F  and  depth-integrated  velocity  field  (U,  V)  over  the  globe  from  the  Poisson  *F-equation.  The 
Poisson  equation- solver  is  similar  to  the  box  model  developed  by  Wunsch.  Thus,  this  method  combines  the  strength  from  both  box 
and  P-vector  models.  The  calculated  depth-integrated  velocity  and  *F-field  agree  well  with  earlier  studies. 
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1.  Introduction 

Winds,  density,  and  friction  determine  volume  trans¬ 
port  in  oceans.  The  wind-driven  volume  transport  has 
been  estimated  using  the  Sverdrup  (1947)  relation  for  the 
interior  ocean,  and  using  the  Stommel  (1948)  and  Munk 
(1950)  linear  frictional  ocean  model  for  the  intensive 
western  boundary  currents.  Leetmaa  and  Bunker  (1978), 
Meyers  (1980),  and  Godfrey  and  Golding  (1981) 
calculated  similarly  for  the  North  Atlantic,  tropical 
Pacific,  and  Indian  Oceans.  Baker  (1982)  examined  the 
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Sverdrup  relation  in  Antarctic  regions.  Godfrey  (1989) 
used  the  Sverdrup  model  with  climatological  annual 
winds  (Hellerman  and  Rosenstein  1983)  to  calculate  the 
mean  depth-integrated  streamfunction  for  the  world  ocean 
under  two  assumptions:  (1)  the  ocean  is  stagnant  below 
some  depth,  and  (2)  all  major  undersea  topographic 
features  such  as  mid-ocean  ridges  lie  below  that  depth. 

The  density-driven  volume  transport  has  been  calcu¬ 
lated  by  several  authors.  The  density  field  directly  deter¬ 
mines  the  geostrophic  velocity  relative  to  the  bottom  flow. 
The  bottom  velocity  ( u-H ,  V-H)  is  usually  calculated  using 
the  /3-spiral  (Stommel  and  Schott,  1977),  Box  (Wunsch, 
1978),  and  P-vector  (Chu,  1995;  Chu  et  al.,  1998a, b;  Chu, 
2000;  Chu  et  al.,  200 la, b;  Chu,  2006)  models. 
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Since  the  geostrophic  balance  is  used,  these  models 
provide  the  solutions  for  the  extra-equatorial  hemi¬ 
spheres,  but  not  for  the  equatorial  region.  An  improved 
inverse  method  is  developed  in  this  study  to  bring 
together  the  two  known  solutions  for  the  extra-equatorial 
hemispheres  across  the  equator  and  to  establish  a  global 
velocity  dataset.  The  rest  of  the  paper  is  outlined  as 
follows.  Section  2  describes  the  basic  theory  of  the 
model.  Sections  3-5  describe  depth-integrated  velocity, 
volume  transport  streamfimction,  and  volume  transport 
vorticity.  Section  6  depicts  the  ‘F-Poisson  equation  and 
its  solver.  Section  7  depicts  the  model  sensitivity.  Section 

8  provides  the  global  circulation  characteristics.  Section 

9  presents  the  conclusions. 

2.  Dynamics 

2.1.  Basic  equations 


The  horizontal  diffusivity  Ah  can  be  estimated  by 
Smargrinsky  parameterization, 

Ah  =  yAxA>i  VV  +  (VV)r|,  (6) 

where 


|VV+(VV)r|  = 


fdv  du\ 
\<9x  dyj 


2 

+ 


Here,  the  nondimensional  parameter  D  varies  from 
0.1  to  0.2  (Mellor,  2003);  we  set  D  =  0.15.  The 
horizontal  grid  in  this  study  is  l°xl°,  i.e.,  (Ax, Ay) 
~  100  km.  Let  the  spatial  variability  of  the  velocity  be 
scaled  by  0.1  m  s-1,  we  have 


|VV+  (VV)Z 


-2  x 


0.1  ms 

105m 


-1 


=  2  x  10” 


(?) 


Let  (x,  y,  z )  be  the  coordinates  with  x-axis  in  the 
zonal  direction  (eastward  positive),  y-axis  in  the 
latitudinal  direction  (northward  positive),  and  z-axis  in 
the  vertical  (upward  positive).  The  unit  vectors  along  the 
three  axes  are  represented  by  (i,  j,  k).  For  the  extra- 
equatorial  region,  the  linear  steady  state  system  with  the 
Boussinesq  approximation  is  given  by 


d2u 

~f{y~v g)  =  Az  ^^2  T"  Ah  V  w, 

(i) 

d2v 

f{u-ug )  =x4z  — +  x4hV2v, 

(2) 

dp 

Tz  =  ~pg ’ 

(3) 

du  dv  dw 
dx~^  dy  +  dz 

(4) 

where  p  is  the  in-situ  density;  /=2L?sin<p,  is  the  Coriolis 
parameter,  Q  the  Earth  rotation  rate,  and  cp  the  latitude. 
V=(w,  v),  is  the  horizontal  velocity;  w  is  the  vertical 
velocity;  V=id/dx+j d/dy,  is  the  horizontal  gradient 
operator;  Vg  =  (ug,  vg),  is  the  geostrophic  velocity 
representing  the  horizontal  pressure  (p )  gradients 

=  _J_dp  =J_d_P  () 

8  f  Pody'  8  f  P,dx'  U 


Substitution  of  (7)  into  (6)  leads  to 


Ah  =  1.5  x  103m2s  1 . 


(8) 


2.2.  Ekman  number 


The  Ekman  number  can  identify  the  relative 
importance  of  the  horizontal  gradient  of  the  Reynolds 
stress  versus  the  Coriolis  force, 


0(14 V2V |)_  A  ,  . 

o{\m  \fF' 

where  L  is  the  characteristic  horizontal  length  scale.  In 
this  study,  the  motion  with  L  larger  than  200  km  is 
considered.  For  extra-equatorial  regions  (north  of  8°N 
and  south  of  8°S), 

\f\>0.2  x  10  4s_I, 


and  the  Ekman  number  is  estimated  by 
1.5  x  103m2s  1 


E<- 


(0.2  x  10'V1)  x  (2  x  105m) 


=  1.875  x  10“3, 


which  shows  that  the  horizontal  gradient  of  the  Reynolds 
stress  can  be  neglected  against  the  Coriolis  force,  i.e., 


Ah  =  0, 


(10a) 


where  p0  is  the  characteristic  value  (1025  kg/m3)  of  the 
sea  water  density.  The  two  coefficients  ( Az ,  Ah )  are  the 
vertical  and  horizontal  eddy  diffusivities. 


for  extra-equatorial  regions. 

For  the  equatorial  regions  especially  near  the  equator, 
|/|  is  very  small.  The  Ekman  number  is  not  a  small 
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Fig.  1.  Bathymetry  of  world  oceans. 


parameter.  The  horizontal  viscous  force,  (AhV2u,AhV2v), 
cannot  be  neglected  against  the  Coriolis  force  in  the 
equatorial  region,  that  is, 


where  z=~H(x ,  y)  represents  the  ocean  bottom,  and 
z=0  refers  to  the  ocean  surface.  Depth-integration  of  (1) 
and  (2)  from  the  ocean  bottom  to  the  ocean  surface  leads 
to 


Ah^0, 

for  the  equatorial  regions. 


(10b) 


JC-rg)=Az^-Af£U_ff+AhvF  {i3) 

—2 AfSJ  u-h  •  •  VH— A\ju-h  V2//, 


3.  Depth-integrated  velocity 

Let  ( U. ,  V)  and  ( Ug ,  Vg)  be  the  depth-integrated 
horizontal  velocity 


dv 

f(U—Us)  =  Az  — 


-A  -\ 

dz  'z=n  zdz'z=~H 


-AbV2V 


(14) 


~2Ab V v-h  •  -VH-Ahv-HS/2H, 


(U,V)=  f°  (u,v)dz, 
J-H 

and  geostrophic  velocity, 


where  ( U-H ,  v-H)  are  velocity  components  at  the  ocean 
(11)  bottom. 

The  turbulent  momentum  flux  at  the  ocean  surface  is 
calculated  by 


(Ug,  Fg)  =  [  (wg,Vg)dz, 
J-H 


(12) 


(du  <9v\  _(zx,Ty) 

\dz'dz)  lz="  _ 


Po 


(15) 
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^=c1 


Fig.  2.  Boundary  conditions  of  'f'for  the  global  ocean. 


where  (tx,ty)  are  the  surface  wind  stress  components. 
The  turbulent  momentum  flux  at  the  ocean  bottom  is 
parameterized  by 


fdu  dv 
Az{  Yz'Tz 


=-H 


=  Cn^uiH  +  v?lh(u-h,v-H), 


(16) 


where  CD  =  0.0025  is  the  drag  coefficient. 

The  thermal  wind  relation  can  be  obtained  from 
vertical  integration  of  the  hydrostatic  balance  Eq.  (3) 
from  the  bottom  (~H)  to  any  depth  (z)  and  then  the  use 
of  the  geostrophic  Eq.  (5) 


Mg  =  U-H  +f-[  -f&z' >  (17) 

f  Pa  J-H  dy 

v‘= (18) 


Substitution  of  (17)  and  (18)  into  the  second  equation 
of  (12)  leads  to 


(t/g,  ve)  =  (t/den  +  Hu-h,  vden  +  Hv-h),  (19) 

where 


( Eden 5  Eden) 


g 

fPo 


d_p 

dx 


dz'dz 


(20) 


is  the  density-driven  transport.  Re-arranging  (13)  and 
(14),  we  have 

AhV2U+fV=fVden+fVb-F  +  AhQu  (21) 

P  0 

-AhV2V  +fU  =fUdea+fUb  +  ^--AhQ2,  (22) 

P  0 

where 

Q\  =  {1S] u-h  •  SIH  +  U-HV2H), 


02  —  (2  V  V-h  •  V//  +  V-tf  V2//), 


and 


Ub  =  (h-  ^  +  u-H,  (23) 

Vh  =  +  -j-  v-H, 

are  the  transport  due  to  the  bottom  currents,  or  simply 
called  the  bottom  transport.  With  the  known  bottom 
velocity  vector  ( u-H ,  V-H),  the  depth-integrated  velocity 
(U,  V)  can  be  determined  from  the  wind,  density,  and 
topographic  data. 
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Fig.  3.  Computed  W- values  for  each  continent/island:  (a)  annual  mean,  (b)  January,  and  (c)  July. 


P.C.  Chu,  C.  Fan  /  Journal  of  Marine  Systems  65  (2007)  376-399 


381 


For  the  extra-equatorial  regions,  the  horizontal  dif¬ 
fusion  can  be  neglected  [see  (8)].  Eqs.  (22)  and  (21) 
become 

u*  =  UAen+Ub  +  fL,  (24) 

JPo 

V*  =  Vden+Vb-)f.  (25) 

J  P  0 

With  the  known  ( u-H,v-H ),  the  depth-integrated  flow 
(£/*,  V*)  may  be  directly  calculated  using  (24)  and  (25). 
However,  the  computed  (£/*,  V*)  field  using  (24)  and 
(25)  is  quite  noisy  and  cannot  not  be  the  final  product. 
Thus,  the  subscript  is  used  to  represent  the  interim 
depth-integrated  velocity  calculated  using  (24)  and  (25). 


is  the  volume  transport  vorticity.  Eq.  (28)  is  called  the 
Poisson  *F-equation. 

5.  Volume  transport  vorticity 

5.1.  Volume  transport  vorticity  equation 


Summation  of  the  differentiation  of  (21)  with  respect 
to  y  and  the  differentiation  of  (22)  with  respect  to  x 
gives  the  volume  transport  vorticity  equation, 


v2n  =  E(K-Fde„-Fb)- 

Ah 


(dQ2_dQi\ 
\  dx  dy  )  ’ 


1 

APo 


(30) 


4.  Volume  transport  streamfunction 


where  P=df/dy,  and  (28)  is  used. 


Integration  of  the  continuity  equation  with  respect  to 
z  from  the  bottom  to  the  surface  yields, 


5.2.  Extra- equatorial  region 


8U  dH  8V  dH 

~7, - F  U-H  — - h  — - h  V-H  —  -W-H  =  0.  (26) 

dx  dx  dy  dy 

With  the  assumption  that  the  water  flows  following 
the  bottom  topography, 


dH 


dH 


W-h  =  U-h  h  V-h  — —  , 

dx  dy 

Eq.  (26)  becomes 

du  dv 
~d7  +  7/-0' 


For  the  extra-equatorial  region,  the  horizontal 
diffusion  can  be  neglected  [see  (8)].  Substitution  of  Ah 
into  0  leads  to 


-7 


d(/Vdea) 


1 

+/L 


dx 

d(fVb) 


dx 


d(/Udenj 


dy 

d{fUb) 


1 


dy  J  /  dx  \p0 


which  leads  to  the  definition  of  the  volume  transport 
streamfunction  *F, 


u  =  - 


dW 


v 


dW 

dx 


(27) 


Subtraction  of  the  differentiation  of  (22)  with  respect 
to  y  from  the  differentiation  of  (21)  with  respect  to  x 
gives 

S/1'F  =  n,  (28) 

where 


1 

~djvdea)_d(fudw)' 

1 

1 

\d{fvb)  _d{fub)-\ 

f 

1 

X 

_ 1 

7 

X 

f  idx  \Po)  dx 
A(d_Qy  d_QA 
+  f\dx  +  dy) 


(29) 


Similarly,  (30)  becomes 

=  (32) 

which  is  the  Sverdrup  relation. 

In  (31),  (C/den,  Eden)  depend  on  p  only;  (rx,Ty)  are 
wind  stress  components;  and  (t/b?  Eb)  are  determined  by 
the  bottom  current  velocity  (u-H,  V-H).  The  P-vector 
inverse  method  (Chu,  1995;  Chu  et  al.,  1998a,b;  Chu, 
2000)  is  used  to  determine  ( U-H ,  v-H)  from  hydrographic 
data  (see  Appendix  A).  In  this  study,  the  climatological 
hydrographic  data  (Levitus  et  al.,  1994)  are  used  to 
compute  (Cden?  Eden)  [see  (20)].  The  climatological 
surface  wind  stress  (rx,Ty)  data  are  obtained  from  the 
Comprehensive  Ocean-Atmosphere  Data  Set  (COADS, 
da  Silva  et  al.,  1994).  The  bottom  topography  is 
obtained  from  the  Navy’s  Digital  Bathymetry  Data 
Base  5-min  (DBDB5)  (Fig.  1).  The  volume  transport 
vorticity  n  is  quite  noisy. 


Latitude  _  Latitude  _  Latitude 
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Fig.  4.  Global  volume  transport  stream  function  ( *P)  computed  from  the  inverse  model:  (a)  annual  mean,  (b)  January,  and  (c)  July. 
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5.3.  Equatorial  region  (between  8°S  and  8°N) 

Let  the  volume  transport  vorticity  n  calculated 
using  (31)  at  8°N  and  8°S  as  the  northern  and  southern 


boundary  values  of  the  vorticity  Eq.  (30).  Here,  the 
forcing  terms  [righthand-side  of  (30)]  are  calculated 
with  the  assumptions  that  (1)  /=/( 8°N)  north  of 
the  equator,  and  /=/( 8°S)  south  of  the  equator,  and 


Longitude 


Fig.  5.  Global  depth-integrated  velocity  ( U ,  V)  vectors  computed  from  the  inverse  model:  (a)  annual  mean,  (b)  January,  and  (c)  July. 
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(2)  ( U. ,  V)  are  calculated  by  (24)  and  (25).  With  the 
given  forcing  terms  and  the  northern  and  southern 
boundary  conditions  and  the  cyclic  eastern  and 
western  boundary  conditions,  the  volume  transport 
vorticity  Eq.  (30)  can  be  solved  in  the  equatorial 
region  between  8°N  and  8°S.  The  computed  El-field  is 
quite  smooth. 

6.  Poisson  -equation 

6.1.  Boundary  conditions 

With  the  computed  global  volume  transport  vorticity 
(IT),  the  Poisson  ‘F-equation  (28)  can  be  solved  when 
the  boundary  conditions  are  given.  No  flow  over  the 
Antarctic  Continent  leads  to  the  southern  boundary 
condition 

W  =  Ci ,  at  the  southern  boundary  y  =  ys.  (33) 

No  horizontal  convergence  of  the  2-dimensional  flow 
( U. ,  V)  at  the  North  Pole  (treated  as  an  island)  leads  to 
the  northern  boundary  condition 

W  =  C2,  at  the  northern  boundary  y  —  yn,  (34) 

where  Cx  and  C2  are  the  constants  to  be  determined. 
The  cyclic  boundary  condition  is  applied  to  the  western 
and  the  eastern  boundaries  (Fig.  2).  We  integrate  8  *F/ 
dy=~U*  with  respect  toy  along  the  western  (or  eastern) 
boundary  from  the  southern  boundary  (*F=0)  to  the 
northern  boundary  to  obtain 

'f'lwestM=-  /V(xwest,/)d/.  (35) 

As 

The  cyclic  boundary  condition  leads  to 

•Pleast  GO  =  ^IwestM-  (36) 

Thus,  the  northern  boundary  condition  is  given  by 

ryn 

¥  =  -  /  U*  (xWGSt,y)dy  =  C2.  (37) 

As 

6.2.  ‘F- Values  at  islands 

Before  solving  the  Poisson  ‘F-equation  (28)  with 
the  boundary  conditions  (33)-(35)  and  (37),  we  need  to 
know  the  ‘F-values  at  all  islands.  These  values  were 


subjectively  set  up  in  some  earlier  studies.  For  example, 
in  calculating  the  geostrophic  transport  in  the  Pacific 
Ocean,  Reid  (1997)  set  up  ‘F-value  to  be  0  for  Antarctic, 
135  Sv  (1  Sv=  106  m3s- ')  for  Australia,  and  130  Sv  for 
America.  In  calculating  the  geostrophic  transport  in  the 
South  Atlantic  Ocean,  Reid  (1989)  set  up  ‘F-value  to  be 
0  for  Antarctic,  132  Sv  for  Africa,  and  130  Sv  for 
America.  Such  a  treatment  subjectively  prescribes 
130  Sv  through  the  Drake  Passage  and  132  Sv  through 
area  between  Africa  and  Antarctica. 

An  objective  method  depicted  in  Appendix-B  is  used 
to  determine  ‘F-values  at  islands.  Fig.  3  shows  the 
distribution  of  ‘F-value  for  each  continent/island 
computed  from  the  annual,  January,  and  July  mean 
hydrographic  and  wind  data.  Taking  the  annual  mean  as 
an  example,  we  have:  0  for  the  American  Continent, 
157.30  Sv  for  Antarctica,  -21.74  Sv  for  Australia, 
-27.17  Sv  for  Madagascar,  and  -21.74  Sv  for  New 
Guinea. 

7.  Model  sensitivity 

With  the  given  values  at  the  boundaries  and  islands, 
we  solve  the  Poisson  ‘F-equation  (28)  with  climatolog¬ 
ical  annual  and  monthly  fJ-fields  and  obtain  annual  and 
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Fig.  6.  Monthly  volume  transport  through  the  Drake  Passage  with  a 
small  seasonal  variation. 
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monthly  global  ^-fields.  After  that,  we  use  (27)  to  re¬ 
compute  the  depth- integrated  velocity  ( U ,  V).  Since 
1°  x  1°  hydrographic  (Levitus  and  Boyer,  1994;  Levitus 
et  ah,  1994)  and  wind  data  (da  Silva  et  ah,  1994)  are 
used  to  compute  ‘F-fields,  small-scale  topographic 
features  such  as  English  Channel,  Taiwan  Strait, 
Gilbralta  Strait,  and  Bering  Strait  cannot  be  resolved 
in  this  study.  Here,  we  present  the  annual  and  monthly 
(January  and  July)  mean  *F-fields  obtained  using  the 
inverse  model.  The  global  *F- field  (Fig.  4)  and  depth- 
integrated  velocity  vector  (£/,  V)  field  (Fig.  5)  agree 
reasonably  well  with  earlier  studies  (e.g.,  Reid,  1989; 
Semtner  and  Chervin,  1992;  Reid,  1994,  1997)  and 
shows  the  capability  of  the  inverse  model  for  determin¬ 


ing  main  characteristics  of  global  circulation  such  as  the 
strong  Antarctic  Circumpolar  Current,  the  well-defined 
subtropical  and  subpolar  gyres,  and  the  equatorial 
current  system.  We  will  discuss  these  features  for  each 
ocean  basin  in  Section  8. 

The  hydrographic  and  wind  data  contain  errors 
(observational  errors).  The  horizontal  diffusion  coeffi¬ 
cient  Ah  is  uncertain.  Sensitivity  study  is  conducted  on 
the  solutions  to  uncertain  Ah  and  the  observational  data 
before  discussing  the  calculated  circulation  character¬ 
istics.  First,  the  model  is  integrated  with  different  values 
of  the  horizontal  diffusivity  Ah.  There  is  almost  no 
difference  among  the  ‘F-fields  with  different  values  of 
diffusivity  Ah  between  1.5  x  103  and  5  x  105  m2  s- 1. 


40°S 


40°S 


Fig.  7.  Inverted  January  F  and  (U,  V)  vector  fields  in  the  Weddell  Sea. 
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Fig.  8.  Inverted  January  •Fand  (U,  V)  vector  fields  in  the  Ross  Sea. 


Second,  suppose  the  observational  data  errors  to  be 
represented  by  a  Gaussian-type  random  variable  (<5%) 
with  a  zero  mean  and  a  standard  deviation  of  cr  whose 
probability  distribution  function  is  given  by 


F{h) 


2  (T2 


(38) 


0.05  N  m-2.  The  model  stability  is  confirmed  from 
almost  no  difference  between  the  ‘F-field  with  random 
errors  in  temperature  and/or  surface  wind  stress  data  and 
the  ‘F-field  with  no  error  added.  Thus,  this  inverse 
model  has  the  capability  to  filter  out  noise  in  the  forcing 
terms  because  of  the  major  mathematical  procedures  of 
the  model  containing  two  integrations  of  the  Poisson 
equation. 


In  this  study,  a  random  number  generator  (FOR¬ 
TRAN  function,  Ranf)  is  used  to  produce  two  sets  of 
random  noises  for  each  grid  point  independently,  with 
mean  value  of  zero  and  standard  deviation  of  cr: 
(a)  three-dimensional  temperature  error  field  with 
standard  deviation  of  0.2  °C  and  (b)  two  dimensional 
surface  wind  stress  error  field  with  standard  deviation  of 


8.  Global  circulation  characteristics 

The  depth-integrated  velocity  vector  field  (£/,  V)  is 
presented  to  illustrate  the  capability  of  this  inverse 
model.  However,  it  is  not  our  intention  to  review  the 
general  circulation  features.  Interested  readers  are 
referred  to  Schmitz  (1996a,b). 
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8.1.  Southern  Ocean 

8.1.1.  Antarctic  Circumpolar  Currents  at  the  Drake 
Passage 

The  computed  monthly  mean  volume  transport 
through  the  Drake  Passage  is  around  156  Sv  with  a  small 
seasonal  variation  (Fig.  6),  which  compares  well  with  the 
estimate  of  1 34  Sv  by  Nowlin  and  Klink  (1986),  although 
observed  Antarctic  Circumpolar  Current  transports  are 
subject  to  uncertainties  of  tens  of  Sverdrups  depending  on 
the  contribution  of  weak  flows  at  depth. 

8.1.2.  Weddell  cyclonic  (clockwise)  gyre 

In  the  Weddell  Sea,  which  probably  contributes  most 
to  bottom  water  formation,  the  water  flows  westward 
under  the  influence  of  the  Coriolis  force  as  it  sinks, 
forming  a  thin  layer  of  extremely  cold  water  above  the 
continental  slope.  It  mixes  with  the  overlying  water, 
which  is  recirculated  with  the  large  cyclonic  eddy  in  the 
central  Weddell  Sea.  The  Weddell  Sea  is  one  of  the  few 
places  in  the  world  ocean  where  deep  and  bottom  water 
masses  are  formed  to  participate  in  the  global  thermo¬ 
haline  circulation.  The  characteristics  of  exported  water 
masses  are  the  result  of  complex  interactions  among 
surface  forcing,  significantly  modified  by  sea  ice,  ocean 
dynamics  at  the  continental  shelf  break  and  slope 
(Foldvik  et  al.,  1985;  Muench  and  Gordon,  1995)  and 
sub-ice  shelf  water  mass  transformation. 

The  most  striking  feature  of  the  inverted  January 
mean  *Fand  ( U. ,  V)  fields  in  the  Weddell  Sea  (Fig.  7)  is 
the  existence  of  the  double-cell  structure  of  the  Weddell 
Gyre  as  suggested  by  the  hydrographic  observations 
(Mosby,  1934;  Deacon,  1979;  Bagriantsev  et  al.,  1989) 
and  the  numerical  simulation  of  a  regional  coupled  ice- 
ocean  model  (Beckmann  et  al.,  1999).  Our  computation 
shows  one  cell  filling  the  western  Weddell  Basin  and  the 
other  trapped  in  a  deeper  basin  northeast  of  Maud  Rise. 
The  west  cell  and  the  east  cell  circulate  4-5  Sv  and 
30  Sv,  respectively.  The  west  cell  is  weaker  than  earlier 
studies  (e.g.,  Beckmann  et  al.,  1999). 

8.1.3.  Ross  cyclonic  (clockwise)  gyre 

The  Ross  Sea  and  adjacent  Southern  Ocean  represent 
important  areas  of  biogenic  production  and  potentially 
large  sources  of  biogenic  material  to  the  water  column 
and  sediments.  Current  meter  moorings  show  that  the 
general  circulation  in  the  Ross  Sea  surface  waters  is 
cyclonic,  with  a  slow  southward  flow  in  the  central  and 
eastern  Ross  Sea  (Pillsbury  and  Jacobs,  1985).  On  the 
basis  of  moored  current  meter  data,  DeMaster  et  al. 
(1992)  pointed  out  the  following  facts:  the  flow  in  the 
southern  Ross  Sea  is  typically  westward.  This  circula¬ 


tion  is  well  developed  in  the  surface  waters  and  extends 
to  depth  as  well.  Current  speeds  at  40  m  above  the 
seabed  are  relatively  low,  on  the  order  of  0. 1  ms-1.  The 
inverted  January  circulation  patterns  (Fig.  8)  shows  that 
the  Ross  Sea  cyclonic  gyre  recirculates  15-30  Sv,  which 
agrees  well  with  Reid  (1997). 

8.2.  Pacific  Basin 

8.2.1.  General  features 

The  annual  mean  northward  transport  across  the 
equator  in  the  west  is  21.7  Sv,  between  the  -21.7  Sv 
isoline  at  the  western  boundary  (northeast  coast  of  New 
Guinea)  to  the  0  Sv  isoline  near  170°W  (Fig.  4a).  This 
current  meanders  and  generates  several  eddies  such  as 
the  Mindanao  Eddy  (cyclonic)  near  southern  Philippines 
and  the  Halmahera  Eddy  (anticyclonic)  near  Indonesia. 
The  northward  current  joins  the  North  Equatorial 
Current  with  30  Sv  transport  (from  0  to  30  Sv  isoline) 
east  of  the  Philippines  (10°-15°N).  Of  these  51.7  Sv  of 
water,  21.7  Sv  are  lost  to  the  Indonesian  seas  directly,  or 
via  the  South  China  Sea  indirectly.  The  remaining  30  Sv 
of  water  continues  northward  to  Japan  and  then 
eastward  with  the  anticyclonic  gyre.  This  subtropical 
gyre  recirculates  20  Sv  between  25°-35°N  (from  30  to 
50  Sv  isoline)  and  makes  50  Sv  of  the  Kuroshio  Current 
east  of  Japan.  Longitudinal  dependence  of  W  along  the 
dateline  and  150°W  (Fig.  9)  in  the  equatorial  region 
show  the  existence  of  the  westward  flowing  South 
Equatorial  Current  (A  *F7  Ay>0)  between  3°S  and  12°S, 
eastward  flowing  Equatorial  Counter  Current  (A  *F/ 
Ay<0)  between  3°S  and  2°N,  and  westward  flowing 
North  Equatorial  Current  (A*F/Ay>0)  between  2°N 
and  12°N. 


Stream  Function  (m2/s)  Stream  Function  (m2/s) 

Fig.  9.  Annual  mean  volume  transport  streamfunction  at  dateline  and 
150°W. 
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There  are  several  low-latitude  cyclonic  gyres,  with 
axes  along  8°N  and  8°S.  Among  them,  an  evident 
cyclonic  gyre  occurs  in  the  north  equatorial  region 
between  180°-120°E,  and  three  smaller  cyclonic  eddies 
(also  called  broken  gyres)  appear  in  the  south  equatorial 
region,  east  of  170°E.  The  north  equatorial  gyre 
appeared  clearly  on  the  study  by  Munk  (1950)  of  the 
wind-driven  circulation  of  the  North  Pacific,  and  was 
identified  by  Reid  (1997)  using  hydrographic  data.  The 
south  equatorial  gyre  identified  in  this  study  (broken 
gyre)  is  different  from  Reid’s  (1997)  and  Tsuchiya’s 


(1968)  description  of  a  complete  gyre  structure.  The 
south  subtropical  anticyclonic  gyre  occurs  east  of  180°E 
between  12°S  and  45°S  and  recirculates  30  Sv  of  water. 

8.2.2.  Pacific-Indian  Ocean  throughflow  region 

The  Indonesian  throughflow  is  the  only  inter-basin 
exchange  of  water  at  low-latitudes  from  the  Pacific  to 
the  Indian  Ocean.  The  calculated  monthly  mean  *F  and 
( U ,  V)  fields  in  the  vicinity  of  Indonesia  (Fig.  10)  shows 
the  volume  transport  and  the  depth-integrated  circula¬ 
tion  pattern  have  weak  seasonal  variations  and  are  quite 
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Fig.  10.  Inverted  annual  mean  F  and  (U,  V)  vector  fields  in  the  vicinity  of  Indonesian  seas. 
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Fig.  11.  Inverted  volume  transport  between  Bali  (8°S,  113.5°E)  and 
northwest  coast  of  Australia  (20° S,  120°E):  (a)  annual  mean  'F- field, 
and  (b)  monthly  variability. 


similar  to  surface  to  intermediate  depth  currents  shown 
in  the  earlier  description  (e.g.,  Fine  et  al.  1994;  Lukas 
et  ah,  1996;  Godfrey,  1996):  water  from  the  Pacific 
Ocean  enters  the  Indonesian  near  the  region  where  the 
New  Guinea  Coast  Current  (10  Sv)  meets  the  Mindanao 
Current  (10-20  Sv),  as  well  through  the  South  Sulu  Sea 
from  the  South  China  Sea  (Fig.  10).  The  New  Guinea 
Coast  Current  transports  10  Sv  water  around  the 
Halmahera  Eddy  into  the  southeastward  flowing  North 
Equatorial  Countercurrent. 

The  North  Equatorial  Current  bifurcates  east  of  the 
Philippines,  with  the  southern  branch  becoming  the 
Mindanao  Current  and  the  northern  branch  becoming 
the  Kuroshio.  Part  of  the  water  flowing  southward  in  the 
Mindanao  Current  retroflects  10-20  Sv  around  the 
Mindanao  Eddy  to  join  the  North  Equatorial  Counter- 
current,  while  the  remainder  enters  the  Celebes  Sea, 
where  some  South  Pacific  water  passes  into  the 
Indonesian  seas  near  Halmahera.  The  majority  of  the 
Indonesian  throughflow  water  exits  to  the  Indian  Ocean 
through  the  Timor  Strait,  with  smaller  transport  through 
the  Savu  Sea  and  Lomok  Strait.  The  Mindanao  Eddy 
(cyclonic)  near  the  southern  Philippines  and  the 
Halmahera  Eddy  (anticyclonic)  near  Indonesia  are 


well  represented  in  the  inverted  W  and  (U,  V)  fields 
(Fig.  10).  No  strong  flow  is  obtained  through  the 
Makassar  Strait  in  the  present  computation  that  is  the 
discrepancy  with  the  earlier  studies. 

The  monthly  volume  transport  between  Bah  (8°S, 
1 13.5°E)  and  northwest  coast  of  Australia  (20°S,  120°E) 
represents  most  of  the  Indonesian  throughflow  (Fig.  11) 
and  shows  a  weak  seasonal  variability  with  a  maximum 
value  of  22.9  Sv  in  December  and  a  minimum  value  of 
20.3  Sv  in  August.  Our  estimation  (Fig.  11)  agrees  well 
qualitatively  with  the  observational  data  (16.2  Sv) 
collected  along  the  same  section  (Bali  to  Australia)  in 
August  1989  by  Fieux  et  al.  (1994)  and  with  the 
numerical  simulated  data  (25.7  Sv  using  the  1.5  reduced 
gravity  model,  15.8  Sv  using  the  nonlinear  six-layer 
model)  reported  by  Morey  et  al.  (1999). 

8.2.3.  Kuroshio  transport  and  its  intrusion  into  the 
South  China  Sea 

Difference  of  the  ‘F-values  between  Japan  and  the 
center  of  the  subtropical  gyre  is  used  as  the  Kuroshio 
volume  transport  (Fig.  1 2).  The  monthly  Kuroshio  volume 
transport  is  very  steady  (very  weak  seasonal  variation) 
with  the  transport  from  57.7  to  62.4  Sv  (Fig.  12).  Our 
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Fig.  12.  Monthly  variation  of  the  Kuroshio  transport. 
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calculation  of  the  Kuroshio  volume  transport  agrees 
with  Schmitz’s  (1996a, b)  estimation  (52.4  Sv). 

The  seasonal  variation  of  the  intrusion  of  the 
Kuroshio  Water  into  the  South  China  Sea  through  the 
Luzon  Strait  has  been  investigated  by  many  authors 
(e.g.,  Shaw,  1989;  Chu  and  Li,  2000).  Shaw  (1989)  used 
the  discriminant  method  to  classify  the  water  mass  T,  S 
characteristics  at  150,  200,  and  250  m,  and  found  that 
the  water  characteristics  of  the  Philippine  Sea  (Kur¬ 
oshio)  were  identifiable  along  the  continental  margin 
south  of  China  from  October  to  January.  The  presence  of 
this  water  indicated  an  intrusion  current  from  the 
Philippine  Sea  into  the  South  China  Sea.  Chu  and  Li 
(2000)  used  the  P-vector  inverse  method  (Chu,  1995, 
2000)  to  determine  the  isopycnal  surface  geostrophic 
velocities  in  the  South  China  Sea.  The  annual  and 
monthly  mean  volume  transports  through  the  Luzon 
Strait  (Fig.  13)  show  Kuroshio  intrusion  all  year  round 
with  a  seasonal  variation  (8-15  Sv).  This  estimate  is 
larger  than  existing  estimations  such  as  2-3  Sv  (Wyrtki, 
1961),  8-10  Sv  (Huang  et  al.,  1994),  2.4-4.4  Sv 
(Metzger  and  Hurlburt,  1996),  and  1.4-13.7  Sv  (Chu 
and  Li,  2000). 

8.2.4.  Australian  Mediterranean  and  South  Australian 
gyre 

The  (*F,  U,  V)  fields  (Fig.  14)  shows  the  following 
features:  the  southward  flowing  East  Australia  Current 
is  the  western  boundary  current  of  the  southern 
hemisphere.  It  is  the  weakest  of  all  boundary  currents, 
carrying  only  about  1 0  Sv.  The  current  first  follows  the 
Australian  coast,  then  separates  from  it  somewhere  near 
34°S  (the  latitude  of  the  northern  end  of  New  Zealand’s 
North  Island).  This  current  recirculates  (10  Sv)  and 
forms  an  anticyclonic  (anticlockwise)  eddy.  The  path  of 
the  current  from  Australia  to  New  Zealand  is  known  as 
the  Tasman  Front,  which  makes  the  boundary  of  the 
warmer  water  of  the  Coral  Sea  and  the  colder  water  of 
the  Tasman  Sea.  In  the  South  Australian  Basin,  an 
anticyclonic  eddy  is  identified  and  recirculates  10  Sv 
with  a  weak  seasonal  variation. 

8.3.  Atlantic  Basin 

8.3.1.  General  features 

Setting  aside  the  part  of  Pacific  west  of  New  Zealand 
and  the  Tonga-Kermadec  Ridge,  which  has  no 
counterpart  in  the  Atlantic,  there  is  some  correspon¬ 
dence  in  the  features  (Figs.  4  and  5).  Each  basin  has  an 
anticyclonic  gyre  in  the  mid-latitudes  of  both  hemi¬ 
spheres,  mostly  west  of  the  major  ridge,  with  two 
cyclonic  gyres  in  between  (low-latitudes). 


The  low-latitude  dual  cyclonic  gyres  have  a  larger 
latitudinal  span  in  the  Atlantic  Basin  (30°S-30°N)  than 
in  the  Pacific  Basin  (10°S-20°N).  The  common  branch 
of  the  dual  cyclonic  gyres  forms  the  eastward  flowing 
equatorial  currents.  For  the  Atlantic  Basin,  there  is  some 
correspondence  between  the  transports  north  and  south 
of  the  equator  aside  from  the  Caribbean  and  Gulf  of 
Mexico,  which  has  no  counterpart  in  the  South  Atlantic. 
The  mid-latitude  anticyclonic  gyres  in  the  North 
Atlantic  and  the  South  Atlantic  have  a  comparable 
strength  (50  Sv).  We  may  also  identify  following  major 
features  from  Figs.  4  and  5:  the  Gulf  Stream,  the 
Labrador  Basin  cyclonic  gyre,  and  the  Brazil-Malvinas 
confluence  zone. 

8.3.2.  Gulf  Stream  volume  transport 

The  Gulf  Stream  volume  transport  can  be  easily 
identified  as  the  difference  of  ‘F-values  at  the  North 
Atlantic  continent  and  at  the  center  of  the  subtropical 
gyre.  The  monthly  mean  Gulf  Stream  transport  is  quite 
steady  with  a  maximum  transport  of  62  Sv  in  October 
and  a  minimum  transport  of  52  Sv  in  March  and  April 
(Fig.  15).  The  calculated  Gulf  Stream  volume  transport 
(57  Sv)  is  too  weak  compared  to  the  value  of  120  Sv 

Annual  Mean  Volume  Transport  Streamfunction  (Sv) 
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Fig.  13.  Monthly  variation  of  volume  transport  through  the  Luzon 
Strait  into  the  South  China  Sea. 
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found  after  detachment  from  Cape  Hatteras  when 
encompassing  the  Southern  Recirculation  gyre  transport 
(Hogg,  1992).  This  is  due  to  the  smoothed  nature  of  the 
climatological  wind  and  hydrographic  data  used. 

8.3.3.  Brazil-Malvinas  confluence 

From  its  littoral  margin  to  the  open  ocean,  the 
western  South  Atlantic  is  marked  by  the  circulation 
patterns  and  exchange  processes  that  are  centrally 
important  to  the  regional  marine  resources  and  local 


economics,  and  equally  important  to  the  global  flux  of 
heat  and  dissolved  substances  (Campos  et  al.,  1995). 
The  depth-integrated  western  boundary  current  (Brazil 
Current)  originates  from  the  South  Equatorial  Current 
(Figs.  4  and  5).  A  major  change  in  the  flow  patterns 
along  the  western  boundary  occurs  in  the  southern 
Brazil  Basin.  Among  other  important  characteristics,  the 
Southwest  Atlantic  is  characterized  by  the  presence  of 
the  Brazil  Current,  a  warm  western  boundary  current 
that,  while  weaker  than  the  Gulf  Stream  in  terms  of  the 


130°E  140°E  150°E  160°E  170°E  180 

Longitude 


Fig.  14.  Inverted  annual  mean  F  and  (U,  V)  vector  fields  in  the  vicinity  of  the  Australian  Mediterranean  sea  and  the  South  Australian  Basin. 
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mass  transport,  is  energetically  comparable  to  its  North 
Atlantic  counterpart,  particularly  in  the  region  of 
confluence  with  the  northward-flowing  Malvinas  Cur¬ 
rent  at  approximately  38°S  (Fig.  16).  The  western  limb 
of  the  recirculation  cell  (anticyclonic)  separates  from  the 
continental  slope  at  about  38°S  upon  its  confluence  with 
the  northward-flowing  Malvinas  Current,  whereupon 
the  bulk  of  the  Malvinas  retroflects  cyclonically 
(clockwise)  back  toward  the  southeast  while  lesser 
portions  continue  northeast  along  the  coast.  On  the 
eastern  side  of  the  cyclonic  trough  is  the  combined 
southeastward  flow  of  Malvinas  and  Brazil  Current 
waters  that  extend  to  45 °S  before  the  subtropical  waters 
turn  east  and  north  to  form  the  pole-ward  limits  of  the 
subtropical  gyre.  The  Malvinas  waters  continue  south  to 
the  southern  margin  of  the  Argentine  Basin  (49°S) 


before  turning  east  with  the  Antarctic  Circumpolar 
Current  regime.  Our  results  are  consistent  with  the 
earlier  studies  (e.g.,  Peterson  and  Whitworth,  1989). 
The  Brazil-Malvinas  confluence  occurs  all  year  round 
with  a  very  weak  seasonal  variability. 

8.4.  Indian  Ocean 

The  depth- integrated  circulation  pattern  in  the  Indian 
Ocean  is  depicted  in  Figs.  4  and  5.  Here,  we  discuss  the 
characteristics  of  the  Agulhas  Current  System  and 
Transport  through  the  Mozambique  Channel  as  exam¬ 
ples  showing  the  model  capability.  The  Agulhas  Current 
is  the  western  boundary  current  in  the  Indian  Ocean.  It 
retroflects,  meanders  and  sheds  discrete  eddies  that 
translate  into  and  across  the  South  Atlantic  Ocean.  This 


Annual  Mean  Volume  Transport  Streamf  unction  (Sv) 


North  Atlantic  Section  Transport 


Monthly 

Fig.  15.  Inverted  monthly  volume  transport  between  the  North  American  east  coast  and  the  center  of  the  subtropical  gyre  representing  the  Gulf  Stream 
transport. 
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is  the  major  contributor  to  the  inter-basin  exchanges  of 
heat  and  salt  between  the  South  Indian  and  the  South 
Atlantic  Oceans  (Gordon,  1985). 

8.4.1.  Agulhas  Current  Retroflection 

The  model  inverts  the  Agulhas  Current  System 
(Fig.  17)  well  with  earlier  depictions  (e.g.,  Gordon, 
1985;  Schmitz,  1996a,b).  The  South  Equatorial  Current 
bifurcates  into  northward  and  southward  branches  at  the 


northwestern  coast  of  Madagascar.  The  southward 
branch  (East  Madagascar  Current)  carries  20-30  Sv 
and  merges  with  the  western  boundary  current  (20- 
30  Sv)  through  the  Mocabique  Channel  to  form  the 
Agulhas  Current  (50  Sv).  This  current  retroflects  at  40°S 
near  the  southern  tip  of  Africa  and  the  return  current 
becomes  the  east  wing  of  a  permanent  eddy  (Agulhas 
Eddy,  anticyclonic)  which  recirculates  10  Sv.  We  also 
see  another  anticyclonic  eddy  (10  Sv)  occurring  in  all 


Fig.  16.  Inverted  annual  mean  F  and  (U,  V)  vector  fields  in  the  southwestern  South  Atlantic  Ocean. 
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Fig.  17.  Inverted  annual  mean  *Fand  (U,  V)  vector  fields  in  the  vicinity  of  the  south  tip  of  Africa  representing  the  Agulhas  Retro  flection. 


monthly  fields  in  the  South  Atlantic  west  of  the  south  tip 
of  Africa. 

8.4.2.  Transport  at  the  Mozambique  Channel 

The  southward  flow  through  the  Mocabique  Channel 
is  a  major  contributor  of  the  Agulhas  Current.  The 
volume  transport  through  this  channel  is  calculated  by 
the  difference  of  ^-values  for  the  African  continent  and 
for  Madagascar.  The  monthly  mean  transport  through 


the  Mozambique  Channel  has  a  weak  seasonal  variation 
with  a  maximum  transport  of  3 1  Sv  in  December  and  a 
minimum  transport  of  23  Sv  in  August  (Fig.  18). 

9.  Conclusions 

(1)  An  inverse  model  is  constructed  to  calculate  the 
depth-integrated  circulation  and  volume  trans¬ 
port  streamfunction  using  surface  wind  and 
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Fig.  18.  Monthly  mean  transport  through  the  Mozambique  Channel. 

hydrographic  data.  This  model  contains  four 
major  components:  (a)  the  P-vector  inverse 
model  to  obtain  the  solutions  for  the  two  extra- 
equatorial  hemispheres  (see  Appendix-A),  (b)  the 
objective  method  to  determine  the  ^-values  at 
individual  islands  (see  Appendix-B),  (c)  the 
Poisson  equation-solver  to  obtain  the  11-values 
over  the  equatorial  region  from  the  volume 
transport  vorticity  Eq.  (30),  and  (d)  the  Poisson 
equation-solver  to  obtain  the  *F  and  depth- 
integrated  velocity  field  ( U. ,  V)  over  the  globe 
from  the  Poisson  ‘F-equation. 

(2)  The  Poisson  equation-solver  is  similar  to  the  box 
inverse  model  developed  by  Wunsch  because  both 
methods  are  based  on  mass  conservation.  The  P- 
vector  method  provides  the  interim  solutions  with 
noise  (appearing  in  the  forcing  terms  in  the 
Poisson  equation)  for  extra-equatorial  region  and 
the  integration  of  the  Poisson  equation  (box  type 
method)  filters  out  the  noise  and  provides  final 
solutions.  Thus,  the  proposed  method  has  the 
strength  from  both  P-vector  and  box  models. 

(3)  This  inverse  model  uses  realistic  topography  and 
has  the  capability  to  filter  out  noise  (in  the  forcing 
terms)  since  two  Poisson  equations  are  integrated. 
The  inverted  volume  transport  is  insensitive  to 


noise  in  wind  and  hydrographic  data.  These  features 
make  it  applicable  for  practical  use. 

(4)  The  inverted  global  and  regional  depth-integrated 
circulation  patterns  agree  well  with  earlier  studies. 
The  monthly  ‘Fand  (U,  V)  fields  provide  realistic 
open  boundary  conditions  for  regional/coastal 
models. 

(5)  An  objective  method  is  developed  to  determine 
‘F-values  at  the  islands  on  the  base  of  Stokes 
circulation  theorem  (see  Appendix-B).  This 
method  contains  two  components:  (a)  an  alge¬ 
braic  equation  for  linking  the  ‘F-value  of  the 
island  to  the  circulation  around  it  and  the  ‘F- 
values  in  the  neighboring  water,  (b)  an  iterative 
algorithm  for  determining  the  ‘F-value  at  the 
island.  Determination  of  ‘F-values  at  islands 
using  this  objective  method  is  not  sensitive  to 
the  noise  level. 
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Appendix  A.  P-vector  inverse  method 

As  pointed  out  by  Wunsch  and  Grant  (1982),  in 
determining  large-scale  circulation  from  hydrographic 
data,  we  can  be  reasonably  confident  of  the  assump¬ 
tions  of  geostrophic  balance,  mass  conservation,  and  no 
major  cross-isopycnal  mixing  (except  for  water  masses 
in  contact  with  the  atmosphere).  Under  these  condi¬ 
tions,  the  density  of  each  fluid  element  would  be  con¬ 
served,  which  mathematically  is  given  by  which 
changes  into 

V  Vp  =  0.  (Al) 

The  conservation  of  potential  vorticity  equation  can  be 
obtained  by  differentiating  (Al)  with  respect  to  z,  using 
geostrophic  and  hydrostatic  balances  (5)  and  (3),  and 
including  the  latitudinal  variation  of  the  Coriolis  parameter, 

V  •  Vq  =  0,  (A2) 

where  q  is  the  potential  vorticity 


It  is  noted  that  neglect  of  relative  vorticity  may 
induce  a  small  but  systematic  error  into  the  estimation  of 
potential  vorticity. 


396 


PC.  Chu,  C.  Fan  /  Journal  of  Marine  Systems  65  (2007)  376-399 


Fig.  Al.  Intersection  of  constant  q  and  p  surfaces  and  the  vector  P. 


When  the  constant  p  surface  intersects  the  constant 
^-surface  (Fig.  Al),  it  is  true  that 


Vp  x  Vg^O. 


A  unit  vector,  called  the  perfect  vector  (or  P-vector) 
by  Chu  (1995),  can  be  defined  by 


V p  x  Vq 
I  Vp  x  Wq\  ' 


(A4) 


Eqs.  (Al)  and  (A2)  show  that  V  is  perpendicular  to 
both  V p  and  V q ,  and  thus,  V  is  parallel  to  P, 

V  =  yP,  (A5) 


where  y  is  a  scalar  and  its  absolute  value  |y|  is  the  speed. 

A  two-step  method  was  proposed  by  Chu  (1995)  (i.e., 
the  P-vector  inverse  method):  (a)  determination  of  the  unit 
vector  P,  and  (b)  determination  of  the  scalar  y  from  the 
thermal  wind  relation.  Applying  the  thermal  wind  relation, 


^  g  r  eP , 

U  =  U-h  +  ~yr~  /  —  dz  , 

>  J-u 


fpo  J-h  dy 


g 


v  =  V-H  +  7 —  /  —  dz' . 


/ 


dp 


fPo  J-H  dx 


(A6) 

(A7) 


at  two  different  levels  zk  and  (or  pk  and  pm),  a  set  of 
algebraic  equations  for  determining  the  parameter  y  are 
obtained 

-yWpW-yWpW  =  A  Ubn,  (A8) 


7«p«-7Mp  (rn)  =  Avfan?  (A9) 

where 


•-/i£ 


»  _  a  . , 


& 


Avkm=~- —  /  —  dz'. 


/Po 


jC 


dp 


dx 


(A  10) 


If  the  determinant 


p(£) 

p(m) 

y 


p(m) 

x 

p(m) 

y 


(All) 


the  algebraic  Eqs.  (A8)  and  (A9)  have  unique  solution 
for  y ^k\m  f  k ), 


y(*)  = 


AUkm  pW 

Avkm  Pf 

pW  p(w) 

pW  p(m) 
y  y 


(A12) 


As  soon  as  y^k)  is  obtained,  the  absolute  velocity  V  is 
easily  computed  using  (A5). 


Appendix  B.  Stokes  circulation  theorem  for 
determining  -values  for  islands 

B.l.  Stokes  circulation  theorem 

Consider  the  (x,  y)  plane  with  uniform  grids 
(Ax, Ay).  Let  the  ocean  basin  be  represented  by 
domains  Q  and  islands  be  represented  by  Qj  (j=  1, 
...,  N)  with  horizontal  boundaries  of  closed  solid-wall 
segments  of  5Qj,j=  1,  ...,  N.  The  domain  Q  may  have 
open  boundaries.  Fig.  A2  shows  a  schematic  illustra¬ 
tion  of  such  a  domain  with  three  open  boundary 
segments,  and  two  islands.  To  determine  the  boundary 
conditions  for  islands,  McWilliams  (1977)  defined  a 
simply  connected  fluid  region  between  an  island  (SQj) 
and  a  clockwise  circuit  in  the  fluid  interior  ( Scoj ).  Let 
Cj  denote  the  closed  area  bounded  by  6Qj  and  dcoj,  and 
(n,  s)  be  the  normal  (positive  outward)  and  tangential 


Fig.  A2.  A  multiply  connected  domain.  The  arrows  indicate  the 
directions  of  integration  along  the  line  integral  paths  defined  in  the 
text. 
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unit  vectors  along  the  boundaries  of  Cy.  The  circulation 
around  the  boundary  of  Cy  is  calculated  using  the 
Stokes  theorem, 

Asnf  'sdl  +  haf  vsxi/ 

k  •  (V  x  V)dxdy,  (Bl) 

where  /  is  the  path  along  the  boundaries  of  Cy,  k  is 
the  unit  vector  in  the  vertical  direction,  V  is  the 
horizontal  gradient  operator.  The  direction  of  closed 
integration  j>  is  defined  as  anticlockwise.  Substitution 
of  (27)  into  the  first  term  in  the  lefthand  side  of  (Bl) 
leads  to 

f>SQjVV*ndI=f5abY*sd! 

(B2) 

k  •  (V  x  V)dxdj/, 

which  links  the  ‘F-values  at  the  island  Qj  to  the 
surrounding  velocity  field.  The  smaller  the  area  of  Cy, 
the  smaller  the  value  of  the  second  term  in  the 
righthand-side  of  (B2),  i.e., 

as  cn(l  (B3) 

where 

r,  =  hcof 


a  circuit  (Scof)  is  called  the  minimum  circuit  along  the 
island  Qj  (Fig.  A3).  Let  (JiJi)  (/=1,..2V+1)  be  the 
anticlockwise  rotating  grid  points  along  Scof  with  (IN+  u 
JN+\)  =  (IU  Ji),  and  let  the  circulation  along  Scof  be 
denoted  by  and  computed  by 


i= l 

+  V(//+1,/m)]  •  [z(//+1 -//) Ax 
+j(JL+l-Ji)Ay\,  (B4) 


which  is  solely  determined  by  the  island  geometry  and 
the  velocity  field  V. 


B.2.  Algebraic  equation  for  W-value  at  island  Qj 
The  lefthand  side  of  (B3)  is  discretized  by 

Fuf'P  ^[nyrl) 

+  'F(Il+l,Jl+1-l)-'F(Il,Jl+  l)-V(Il+uJl+l  +  1)] 

i=  l  ^ 

+  l,//+1)-^/-l^/)-^/+i-l^m)]- 


Thus,  selection  of  Scoj  with  a  minimum  Cy  becomes  a 
key  issue  in  determining  the  streamfimction  W\ q.  Such 


Since  the  grid  points  on  the  island  Qj  are  always  on 
the  left  side  of  the  anticlockwise  circulation  fj  (Fig.  A3), 


Fig.  A3.  Grid  points  of  the  minimum  circuit  along  the  island  Qj. 
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half  grid  points  of  (A4)  are  in  the  island  and  half  in  the 
water.  Thus,  Eq.  (B4)  can  be  written  by 

r}=AT\ Qj  +  rf 

where  i}(w)  is  the  circulation  in  the  water  and 

,  _  ^  f\Ji+x~Ji\Sy  , 

h  V  2Ax  +  2Ay  )■ 

The  volume  transport  streamfunction  at  island  Qj  is 
computed  by 

=  (B5) 


B.3.  Iteration  process 

Eq.  (B5)  cannot  be  directly  used  to  compute  *F| Q. 
even  if  the  vertically  integrated  velocity  (U,  V)  is  given. 
This  is  because  that  the  ‘F-values  at  surrounding  water  is 
still  undetermined.  Thus,  we  need  an  iterative  process  to 
determine  *F| q.  from  a  first  guess  value. 

Suppose  all  the  islands  Qj  (j= 2, ...,  N)  in  Fig.  A2  to 
be  removed.  With  the  given  boundary  conditions  at 
we  solve  the  ‘F-Poisson  equation  (28)  and  obtain  the 
solution  *F*(x,  y).  Average  of  ‘F*  over  Qj  leads  to  the 
first  guess  ‘F-values  at  islands  Qj  (j= 2, ...,  TV), 


?V.{0)  =JJq  >P*(x,y)dxdy. 


Let  ‘F-values  and  the  circulation  be  given  at  the 
m- th  iteration  such  that 


Via  (m)  = 


r, 


(B6) 


where  the  minimum  circuit  circulation  at  the  m- th 
iteration,  rf(m ),  might  not  be  the  same  as  /}.  We 
update  *F\q  using 

Subtraction  of  (B6)  from  (B7)  leads  to 

^\0j(m+l)  =  V\Qj(m)+  j  , 

which  indicates  the  iteration  process:  (a)  solving  the  ‘F- 
Poisson  equation  (28)  with  *F| aftri)  to  obtain  solutions 


(B7) 


(B8) 


and  in  turn  to  get  2}*(w),  (b)  replacing  the  ‘F-values  at 
islands  using  (B8).  The  iteration  process  repeats  until 
reaching  a  certain  criterion 


|  dr* 

Trj 


-<e, 


(B9) 


where 


\r\= 


\ 


(BIO) 


and  e  is  a  small  positive  number  (user  input),  which  is 
set  to  be  1(T6  in  this  study.  As  soon  as  the  inequality 
(B9)  is  satisfied,  the  iteration  stops  and  the  final  set  of 
{  *F|^,  /=  1 ,  2,  ...,  N}  become  the  optimal  ‘F-values  for 
islands. 
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